Edge Current due to Majorana Fermions in Superfluid 3 He A- and B-Phases 
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We propose a method utilizing edge current to observe Majorana fermions in the surface Andreev 
bound state for the superfluid 3 He A- and B-phases. The proposal is based on self-consistent analytic 
solutions of quasi-classical Green's function with an edge. The local density of states and edge mass 
current in the A-phase or edge spin current in the B-phase can be obtained from these solutions. The 
edge current carried by the Majorana fermions is partially cancelled by quasiparticles (QPs) in the 
continuum state outside the superfluid gap. QPs contributing to the edge current in the continuum 
state are distributed in energy even away from the superfluid gap. The effect of Majorana fermions 
emerges in the depletion of the edge current by temperature within a low-temperature range. The 
observations that the reduction in the mass current is changed by T 2 -power in the A-phase and 
the reduction in the spin current is changed by T 3 -power in the B-phase establish the existence 
of Majorana fermions. We also point out another possibility for observing Majorana fermions by 
controlling surface roughness. 

PACS numbers: 67.30.hp 



I. INTRODUCTION 

The superfluid 3 He is firmly established as an 
anisotropic p-wave superfluid [1| . At an edge of the su- 
perfluid 3 He, the surface Andreev bound state emerges 
by the pairing symmetry due to the bulk-surface corre- 
spondence. Recently, the surface Andreev bound state 
has become a focus of concern from topological aspects. 
The superfluid gap of the superfluid 3 He among topo- 
logical superfluids is closed at the interface of a topolog- 
ically trivial vacuum by a topological phase transition. 
This is a topological aspect of the surface Andreev bound 
state 0]. The topological features arc quite different be- 
tween the superfluid 3 He A- and B-phases because the A- 
(B-)phase is a time reversal symmetry broken (unbroken) 
state. The A-phase is a chiral superfluid with the sponta- 
neous edge mass current while B-phase is a helical super- 
fluid with the spontaneous edge spin current. Quasipar- 
ticles (QPs) bound in the surface Andreev bound state 
behave as Majorana fermions owing to the particle-hole 
symmetry [3, [4|. The Majorana nature in the A- and 
B-phases is also distinctive. Majorana fermions have a 
linear dispersion relation forming a "Majorana valley" 
in the A-phase 0, @] and a "Majorana cone" in the B- 
phase 0, Q . This Majorana cone has been observed by 
recent experiments [3, [lQl ■ However, there has been no 
firm evidence of the edge current that accompanies Ma- 
jorana fermions. 

For other candidates of the topological supercon- 
ductor, e.g Sr 2 Ru0 4 [H 0, UPt 3 [H, [lj, and 
Cu K Bi 2 Se 3 [HljGil: the precise pairing symmetry has not 
yet been identified and has been under intense discussion. 
Thus, there has been no concrete topological supercon- 
ductors established as the superfluid 3 He in the chiral 
A-phase and helical B-phase thus far. 

In the two-dimensional (2D) Fermi surface model of 
the A-phase, namely, the fully gapped chiral p-wave state 
as the proposed pairing symmetry for Sr2Ru04 [111 ], the 



edge mass current is carried both by Majorana QPs 
within the superfluid gap and by the QPs in the con- 
tinuum state outside the gap [IT], 03 ■ These two kinds 
of QPs have different contributions to the edge mass cur- 
rent and angular momentum. Stone and Roy [l7( first 
proved that, in the 2D A-phase under a uniform pair 
potential, the magnitude of the angular momentum by 
the edge mass current is Nh/2 at a zero temperature (N 
is the total number of 3 He atoms in the whole system). 
The A-phase with the three-dimensional (3D) Fermi sur- 
face under the pair potential reflected on an edge has the 
same magnitude of this angular momentum and the same 
contributions from these two kinds of QPs [19| . 

It is also expected that Majorana QPs and QPs in the 
continuum state have different contributions to the edge 
spin current in the B-phase. In this study, we investigate 
the roles of these two kinds of QPs in the edge spin cur- 
rent with the density of states (DOS) for these QPs. It 
turns out that the edge spin current from Majorana QPs 
is partially cancelled by that from QPs in the continuum 
state similarly to the A-phase [l9| . Finally, the total spin 
current at a zero temperature is J s = — (K/27r)(nft/6), 
where k is the quantum of circulation and n is the den- 
sity of 3 He atoms. The numerical coefficient of nh for 
the spin current is 2/3 of that for the total mass current 
J = —nh/4 in the A-phase because QPs contributing to 
the spin current are in two of three spin states of QPs. 
The temperature dependence of the spin current will also 
be shown. The depletion of the spin current by temper- 
ature within a low-temperature range is proportional to 
T 3 because the spin current from Majorana QPs has a 
quadratic energy spectrum. 

We base our arguments on the quasi-classical theory, 
which is valid for £ ^$> kp 1 , well satisfied for the superfluid 
3 He (coherent length £ ~ 10-100 nm and Fermi wave 
number kp ~ 0.1 nm). We find analytic solutions for 
the pair potential when the system has an edge. These 
analytic solutions give useful and transparent informa- 



tion on various physical properties at the edge. Moreover, 
these solutions are self-consistent solutions in the infinite 
cutoff energy limit, namely, the "weak- coupling limit". 
Thus, the present comparative studies on both the A- 
and B-phases based on the same theoretical framework 
may help us understand the nature of Majorana fermions. 
This paper is arranged as follows: In i jlll we formu- 
late the quasi-classical theory on the basis of the quasi- 
classical Green's function, which gives quantitative infor- 
mation on QPs. We discuss the edge mass current and 
angular momentum in the A-phase in ^IIII in compari- 
son with the 2D Fermi surface model. We also touch on 
the so-called intrinsic angular momentum problem. In 
^IV[ we investigate the role of Majorana QPs in the edge 
spin current with the DOS in the B-phase. We also show 
the temperature dependence of the spin current. We de- 
vote the final section to the summary and discussion. We 
demonstrate the self-consistency of our solutions in the 
weak-coupling limit in Appendix. 



II. QUASI-CLASSICAL THEORY 

Microscopic information on the edge state is contained 
in the quasi-classical Green's function g(r,k,u n ). The 
quasi-classical Green's function is calculated using the 
Eilcnberger equation [2(| as 

-ihv(k) ■ Vg{r, fc, w„) = iuj n a z - A(r, fe), g(r, k, u n ) 

(1) 

In this paper, the "ordinary hat" and "wide hat" indi- 
cate the 2x2 matrix in spin space and the 4x4 ma- 
trix in particle-hole and spin spaces, respectively. The 
quasi-classical Green's function is described in particle- 
hole space by 



g{r,k,u n ) = -i-n 



g(r,k,w n ) if(r,k,uj n ) 
-if(r,k,u n ) -g(r,k,(jj n ) 



(2) 



with the center-of-mass coordinate of a Cooper pair r, 
the direction of the relative momentum of a Cooper pair 
k, and the Matsubara frequency uj n = (2n + Vf-nksT 
with n £ Z. The quasi-classical Green's function satisfies 
a normalization condition g 2 = — 7r 2 l. The pair potential 
is described in particle-hole space by 



A(r,fc) 



A(r,fc) 

-A+(r,fc) 



(3) 



and that in spin space for the spin-triplet state is defined 
by the d-vector as 



A(r, fe) = (icrcTy) ■ d(r, fe), 



(4) 



with the Pauli matrix a. The Fermi velocity in the Eilcn- 
berger equation is given as v(k) = vpk for the 3D Fermi 
sphere. 



We solve eq. CD b y the Riccati method 12lH23j. We 
introduce the Riccati amplitude a = (1 + g)^ 1 f and b = 
(I+5) -1 / related to particle- and hole-like projections of 
off-diagonal propagators, respectively. Equation (jTJ) can 
be rewritten as Riccati equations: 

hv(k) ■ Va(r,fc,w n ) 

= A(r,fe) -a(r,fc,w„)A t (r, k)a(r,k,u n ) 

- 2u n a(r,k,u n ), 
— hv(k) ■ V6(r, fc, ui n ) 

= A 1 (r, fc) - b(r, fe, w„)A(r, fc)S(r, fc, w„) 

-2ojJ)(r,k,uj n ). (5) 

These equations are solved by integration toward fc for 
a(r, k,uj n ) and toward — fc for b(r, k,uj n ). By the Riccati 
amplitude, the quasi-classical Green's function is given 
as 



g = -iTT 



(l + ab)- 1 







(l + ba)- 1 ; V —2*6 -(1 -6a) 



1 — ab 2ia 



(6) 



By using the quasi-classical Green's function, the 
temperature-dependent mass and spin currents are cal- 
culated using 

j(r,T)=mN 7rk B Tj2(v(k) Im [g {r, fc, w„)]> fc) (7) 
j»{r,T) = ~N irk B Tj^(v(k) Im^(r,fc ) w„)])i ) (8) 



respectively, where N is the DOS in the normal state, m 
is the mass of the 3 He atom, (■ ■ -)k indicates the Fermi 
surface average, and g^ is the component of the quasi- 
classical Green's function g in spin space, namely, 



9 



.90 + 9z 9x - ig v 

fjx + ig y go - g z 



The mass and spin currents and local density of states 
(LDOS) for the energy E are given by 



j(r,E) = (j(r,k,E)) k 

= mN (v(k) Re[g (r,k,u) n )\ iUn ^E+in]) k ■ 

j»(r 1 E) = (j»(r,k,E)) k 



"(9) 



= -N (v(k) Rfi[flf M (r,fe,w n )|i Wn _ f B+i t j]) fc , 

(10) 
N(r,E) = (N(r,k,E)) k 

= N (Rc[g (r,k,uj n )\ iLOn ^E+i V }) k , (H) 

respectively, where r\ is a positive infinitcsimally small 
constant. 



III. A-PHASE 

In this section, we consider the superfhhd 3 He A-phase 
in a slab with a small thickness along the z-dircction. In 
the sufficiently thin slab, the /-vector pointing to the di- 
rection of the angular momentum of a Cooper pair is 
aligned toward the z-direction and the d-vector is also 
aligned toward the z-direction by the dipole interac- 
tion QJ. This system is realized using a slab of sub- 
orn thickness [6J, which is already realized experimen- 
tally [241 ] . The uniform pair potential in the system is 
described by the d- vector as 



d(k) = A A (k x +ik y )z, 



(12) 



where z is a unit vector and A ,4 is the amplitude of the 
supcrfluid gap in the A-phase. 

Here, we discuss a side edge of the slab at x = which 
is filled with the superfiuid 3 He at x > 0. Assuming that 
the side edge is specular, only the fc^-component of the 
d-vector is suppressed at the edge and the d-vector must 
recover to the chiral state in the bulk as eq. (|T2l far from 
the edge. Accordingly, we solve the Riccati equations in 
eq. ([5]) under the pair potential with 



d z (x, k) = A j 



k x tanh 



u 



i %K% 



(13) 



where the coherent length is defined by £a = Hvf/Aa- 
This d-vector form embodies the fact that the k x - and 
fc y -components are suppressed and intact at the edge, 
respectively. Since the pair potential in spin space is 



A = 



4 

d z 



we can obtain the Riccati amplitude in spin space as 



o\ f _ / b 
a ' U 



Thus, we can reduce the matrix Riccati equations in 
eq. ([5]) to the scalar Riccati equations 



d 

frvFk x —a =d z 
ox 





d* z a 2 — 2u) n a, 



-hv F k x —b =d* 



d z b 2 - 2uj n b. 
The solution of the Riccati amplitude is 

a(x, k,u n ) = 
b{x, k,u) n ) = 



(14) 



u„ 


-yfi% 


+ A A 


sin 2 


9 + d z (x 


k) 


u„ 


+ v5 


+ A A 


sin 2 


9- 


-d* z (x 
-d* z (x 


k) 


u„ 


-\fi% 


+ A\ 


sin 


9- 


k) 



, (15) 



/Ljl + A 2 A sin z e + d z (x,k) 

cos cj) and k y = sin 9 sin <f> i 
spherical coordinates on the unit Fermi surface. This 



where we use k x — sin 9 cos cj> and k y = sin 9 sin <f> in the 



solution formally satisfies the boundary condition [2f 
namely, 



o(— co, k,ui n ) 
b(+oo, k,u> n ) 



d z (-oo,k) 



w„ + y/u% + |d 2 (-oo,fc)| 2 ' 
d* z (+oo 7 k) 

W„ + y/u}2 + \d z ( + 00,k)\ 2 ' 



(16) 



Moreover, this is a self-consistent solution in the weak- 
coupling limit as demonstrated in Appendix. The quasi- 
classical Green's function is obtained using the relation 
g = (1 + dfe) _1 (l — ab). Finally, the spin component of 
quasi-classical Green's function go = (1 — ab)/(l + ab) is 
only finite in g as 



g (x,k,uj n ) 



A\ sin 2 9 



A\ sin 2 9 cos 2 4> 
2(cj„ + i A a sin 9 sin < 



scch 



(17) 



A. Local density of states 

From the quasi-classical Green's function in eq. (|17l) . 
we can calculate 0-angle-resolved LDOS using eq. (fTTj) as 



N(x, 9,E)= j —N(x, k,E) = — scch 2 



for the Majorana bound state \E\ < A A sin9 and 



(18) 



N(x,9,E) = N a 



\E\ 



E 2 - A\ sin 2 1 




\E\ 



E 2 - A\ sin 2 9 



1 sech 



(19) 



for the continuum state \E\ > A^sinS. Thus, QPs feel 
the pair potential A^ sin 9 where 9 is the polar angle from 
the point node situated at the pole of the Fermi sphere. 
The 6*-angle-resolved LDOS at the edge x = for 9 = tt/2 
is shown in Fig.fJJa). The LDOS from the Majorana zero 
energy mode has a constant and finite value, No/2. A 
similar value is also obtained by numerical calculation in 
a finite temperature [a, |f|. The ^-angle-resolved LDOS 
from Majorana QPs is independent of the polar angle 9. 
This indicates that the dispersion of Majorana QPs forms 
a "Majorana valley" [a, III- The LDOS in the continuum 
state is also deformed to conserve the DOS against the 
appearance of Majorana QPs. 



B. Edge mass current and angular momentum 

We can also calculate S-anglc-resolved mass current 
along the edge using eq. © as 



j y (x,e,E)=J^jy(x,k,E) 



mv F N E ,2 f x . 
^A~ A S6Ch [U ] (20) 



for the bound state \E\ < Aa sin# and as 



j v (x,9,E) = 



mvpNo E 



E 2 - A\ sin 2 1 



E 2 



\E\ 



E 2 - A\ sin 2 i 



sech ( — - 



(21) 



for the continuum state \E\ > A^sin^. The #-angle- 
resolved mass current at the edge for 6 = it/2 is shown 
in Fig. QJb) . The energy spectrum of the mass current 
from Majorana QPs is linear. The edge mass current 
from the continuum state is due to the deformation of 
the LDOS, which decreases away from A^sinfl. The 
asymptotic behavior of the edge mass current in eq. (|21[) 
at x = is estimated as 



jy(x = 0,d,E)n ^— sin ( 



A a sin ( 
E~ 



(22) 



for \E\ 3> A^sin^, implying that the contribution of 
the QPs to the edge mass current decreases as ~ E~ 3 . 
This power functional behavior of the decrease means 
that QPs contributing to the edge mass current are not 
confined in only the vicinity of the Fermi surface. 

Since QPs fill the energy state up to the Fermi energy 
at a zero temperature, the mass current along the edge 
from the Majorana bound state is obtained as 



3f\x) = 



dEj y (x,k,E) 

A a sin 9 I fc 

6 S6Ch I U) ' (23) 



and that from the continuum state is obtained as 



fy° nt (x) = 



-Aa sin 6 



mvpNoAj 
12 



dEj y (x,k,E) 



sech ( — - 



6 



(24) 



The total mass current along the edge from the Majorana 
bound state is 



./, 



MJ 



dxjf\x) 



nh 



(25) 



s€ 


(a) 1 










Q5 1 


— J\ 


^ 


^ 0.5 
n 


j ! ' 




|0.5 

c£ 

t= 
4.0.5 

o 
•»1 


(b) 


r 



-A A A A 
E 



FIG. 1. (Color online) Energy profiles of 6-angle-resolved 
LDOS (a) and mass current along the edge (b) at x — for 
9 = 7r/2 in the A-phase. At a zero temperature, QPs fill the 
colored (shaded) states in (a). The mass currents from the 
bound and continuum states are derived by integrating the 
blue (light gray) and pink (gray) regions in (b), respectively. 



and that from the continuum state is 



dxf y 



**(x) - -f, 



(26) 



where the density of 3 He atoms n emerges from the nor- 
mal DOS N = {Z/mv 2 F )n. The same result with eq. (gSJ) 
is obtained in connection with the chiral superconductor 
Sr2Ru04 [26|. Since these currents flow oppositely, the 
total mass current induced by the edge state is 



J y J : 



MJ 



Tcont _ n ' i 

J y ~ 4 • 



(27) 



In a disk with a large radius E>^, since the mass cur- 
rent can be regarded as localized at the edge, the angular 
momentum from each state is calculated as 



L 



MJ 



Nh, L\ 



Nh 



(28) 



Finally, the total angular momentum simply becomes 



-MJ 



+ L<; 



Nh 



(29) 



The cancelled angular momentum Nh/2 may be related 
to our finding that the DOS of Majorana QPs is No/2. 
The total angular momentum due to the edge mass cur- 
rent coincides with that in the 2D A-phase under a uni- 
form pair potential [13, LUt- Interestingly, half of the 
angular momentum from the Majorana bound state is 
canceled by that from the continuum state. 



C. Temperature dependence 

The temperature dependence of the angular momen- 
tum L Z (T) by the edge mass current is calculated using 
cq. as 



L Z (T) -- 

*E 



4 A A 



Bui 



Ai 



A.. 



Ai 



sm 



\A4 + 



Ai 



A, 



(30) 



This temperature dependence is shown in Fig. [2{a) by 
open circles. The component of the supcrfluid density 
tensor parallel (perpendicular) to the direction of the 
point nodes p°,, (p® ± ) [27| is also depicted by two lines. 
The angular momentum L Z (T) (open circles) has the 
same temperature dependence of p®»(T) (solid line), as 
pointed out by Kita |28|. This complete correspondence 
of the temperature dependence, however, may be acci- 
dental because the low-tcmpcraturc depletion of the an- 
gular momentum can be explained by only the Majorana 
QPs, as mentioned below. Moreover, in the 2D Fermi sur- 
face model, the temperature dependence of the angular 
momentum has no connection with that of the supcrfluid 
density [see Fig.[2jb)]. 

According to eq. (|30|) . we can derive the low- 
temperature behavior of the angular momentum as 



L,(T) = f 



irksT 







irksT 



(31) 



The low-temperature depletion of the angular momen- 
tum contains both contributions of the excitations of Ma- 
jorana QPs and point nodes. Majorana QPs contribute 
to the depletion as that in T 2 -power because j y (x, 9, E) 
in the Majorana bound state has a linear energy depen- 
dence [Fig. QJb)]. The contribution of the excitations at 
point nodes is derived from the low-energy behavior of 
the energy spectrum of the mass current in the contin- 
uum state. The energy spectrum at x = is obtained 
using eq. (|2"Tj) as 



jy(x = 0,E) 



m0(ffljy(x = O,6,E) 




for \E\ < A A , where 9 = sin -1 (|.E|/A, 
energy behavior is estimated as 



j„(x = 0,£)w- 



mvpNQ 
12 



E 

a7 



This low- 



(33) 



for \E\ <C Aa- Thus, the excitations at point nodes 
contribute to the angular momentum as the fourth order 
of temperature (~ T 4 ). Thus, the observation of the 
depletion of the angular momentum as that in T 2 -power 
could establish the existence of Majorana QPs. 

In the 2D Fermi surface model, where point nodes are 
absent, LDOS N(x, E) is obtained by the substitution of 
A 2 d for A^sin^ on the right-hand side of cqs. (fT5|) and 
(|i"9|) , where A 2 d is the amplitude of the superfiuid gap 
in the 2D chiral p-wave state. The values of the total 
mass current and angular momentum are also given by 
eqs. ([2"5 ]) -([2"7 ]l and eqs. ([28) and (|2"5|). respectively. How- 
ever, the temperature dependence of the angular momen- 
tum is different from eq. ()30[) . that is, it is given by 



L?(T) 



E 



irksT 



i "" A 2D 




V"l + A 2D < 


oKI 


A 2D A 2D ^ n + Ai D 


"A 2D 



(34) 



This temperature dependence is also shown in Fig. ^h) 
with the superfiuid density in the 2D chiral p-wave state 
P°s2D- These behaviors are clearly different as pointed 
out by Sauls [18[ because the low-temperature depiction 
of the angular momentum is due to the excitations of 
Majorana QPs. The low-temperature behavior of the 
angular momentum is 



2D 



(D - f 



_2 ( irk B T 
3 I A 2D 







irkfiT 



±2D 



(35) 



The energy spectrum of the mass current j y (x, E) in the 
Majorana bound state is obtained by the substitution of 
A 2 d for A a on the right-hand side of eq. (|20|) . The dif- 
ference between the coefficients on the second order of 
temperature in eqs. (f3~Tj) and (|3"5|) comes from the dif- 
ference in the normal DOS, where A^o = (2/m,Vp)n for 
the 2D Fermi surface is 2/3 of that for the 3D Fermi 
sphere. Although the coefficients are different, the tem- 
perature dependences of the angular momentum are sim- 
ilar in these models owing to A a > A 2 £>. 



D. Intrinsic angular momentum 

The angular momentum by the edge mass current at 
a zero temperature is L = Nh/2. This magnitude corre- 
sponds to one of the predictions of the intrinsic angular 
momentum [29l - l3l| . The magnitude of the intrinsic an- 
gular momentum in this prediction is expected naively 
when all Cooper pairs carry one unit of angular momen- 
tum. However, the edge mass current is carried by a 
portion of QPs because the contribution of QPs to the 
edge mass current decreases as E~ 3 away from the su- 
pcrfluid gap. QPs carrying the edge mass current are 



1.0 
0.8 
0.6 

0.4 

0.2 



0.8 
0.6 

0.4 

0.2 



(a) 



L z (T)/(Nh/2? 

■P°s\i(T)/ P 
■P° S ±(T)/P 




Ll D (T)/(Nh/2) 
■P°s2d(T)/ P 



0.4 0.6 

T/T c 



FIG. 2. (Color online) Temperature dependence of angular 
momentum (open circles) in the A-phase (a) and 2D chiral p- 
wave state (b) with the superfluid density. The components of 
the superfluid density tensor p° g « (solid line) and p® ± (dotted 
line) are shown in (a) and the superfluid density p® 2 D (solid 
line) is shown in (b). L Z (T) and p°,i(T) in (a) completely 
coincide. 



dipolc interaction when 9d = 9l = cos _1 (— 1/4) and the 
n- vector is perpendicular to a surface [l| . 

In this section, we consider the superfluid 3 Hc B-phasc 
filled in z > with a specular edge at z = 0. Only 
the /c 2 -component of the d-vector is suppressed near the 
edge within the coherent length and the d-vector must 
recover to the bulk form in eq. (|36|) far from the edge. In 
this system, the n-vector points to the z-dircction per- 
pendicular to the edge, and the angle 9d remains intact 
at the bulk value 9l because the angle can vary not on 
the order of the coherent length but over a larger length, 
namely, the dipole coherent length [l[. Accordingly, the 
pair potential with a specular edge at z — is 



A(z,k) = A B 



/-sm9e~' l ^ +e ^ cos(9tanh(z/£ B ) 
\cos9tanh(z/£ B ) sm9e l ^+ ^ 



(37) 



where we use k x = sin 9 cose/), k y = sin 9 s'mcj), and k z = 
cos 8 and the coherent length is defined by £, B = Ivuf/ 'A^. 
The Eilenbcrgcr equation eq. ([T]) can be solved sim- 
ilarly to the A-phase after the unitary transformation 
using the unitary matrix [32j 



M 



M 0W l/« u* 
M*J ' y/2 \ u ~u* ' 



(38) 



where u 2 = —ie l ^ +0L \ The unitary-transformed pair 
potential is 



A' = MAM f 



A' 

-A't o 



(39) 



distributed neither in the narrow energy shell around the 
Fermi level nor up to the bottom of the Fermi sea. The 
total number of 3 He atoms N in the angular momen- 
tum by the edge mass current comes from the normal 
DOS A^o = (3/mv F )n. Moreover, the temperature de- 
pendence of the angular momentum by the edge mass 
current has no connection with that of the superfluid den- 
sity, which is related to the number of Cooper pairs. The 
low-temperature depletion of this angular momentum as 
that in T 2 -power is due to the excitations of Majorana 
QPs in the edge state. 



IV. B-PHASE 



where 



A'=MAM T -f d ' (z,fl) ° 

-d'*{z, 



with 



d'(z,6) = A B 



cos#tanh | — | 



(40) 



(41) 



Finally, the unitary-transformed Eilenbcrgcr equation is 



i) 



ihvp cos 6*— —g'(z,9,uj n ) 
oz 



= iuj n d z -A'(z,9),g'(z,9,uj n ) , (42) 



The uniform pair potential in the superfluid 3 He B- 
phase is described by the d-vector as 

d{k) = A B R(n, 9 d )(k x x + k y y + k z z), (36) 

where x, y, and z are unit vectors, A B is the amplitude 
of the superfluid gap in the B-phase, and R(n, 6d) is a 
rotation matrix with a rotation axis n and a rotation 
angle 9d about n. The rotation matrix gives the relative 
angle between k and d. The spin state is stable by the 



with g' = M1)M\ where M and a z are commutable. 

Since the unitary-transformed pair potential in spin 
space is described by eq. (|40]) , we can obtain the Riccati 
amplitude in spin space: 




-a' 



, y = 



b 1 
-b" 



Owing to the analogy of the A-phase in i]IIIl the solution 
of the Riccati amplitude is obtained by the substitution 



of Ab and d! for A^ sin 9 and d z , respectively, as 



a'(z,6,u n ) = 
b'(z,9,uj n ) = 



Un 


-V 


-£ + A| 


+ d'(z, 


0) 






-d'*(2 
-d"(z 


,0) 
0) 



w„ + v /o;2+A|+d'(z, 



(43) 



The quasi-classical Green's function is obtained by the 
unitary transformation of this solution as g — M^g'M. 
The obtained quasi-classical Green's function is self- 
consistent in the weak-coupling limit similarly to the A- 
phase case, as discussed in Appendix. The spin com- 
ponents of quasi-classical Green's function g — M'g'M 

I 



,9o = Re 



1 - a'b' 



9x 



1 + a'b' 
g y = —i cos(<p + 6l) Im 

Since 

1 - a'b' 1 



1 + a'b' y/uji + A 



B 

A 2 B cos 2 1 



= ism(cf) + 9l) Im 

I -a'b'] 

, , ,,, , 9z = 0. 
1 + a'b' 



1 - a'b' 



1 + a'b' 



sech 



(44) 



2(u> n + iA B sin9) ~ \^b 

by the analogy of the A- phase in eq. (fTTf , we obtain 



(45) 



g (z,k,u n ) 



V^f+Al 



A| cos 2 9 



g x (z,k,u n ) = 



sin((j)+9 L ) A|jCos 2 6> 
V^ + Al" 



g y (z,k,u n ) 
g z (z,k,u n ) = 0. 



4 
cos(4>+9 L ) A 2 , cos 2 1 



\A4+A 



1 1 \ , 2 /z 

nt : — ;H ^ : — « sech — - 

L0 n + iA B sm9 Lu n -iA B sm9J \£ B 



1 



1 



cjn + iAssinfl w n — iAssinf? 
1 1 



uj n +iA B sin9 w n — iAssmfl 



sech I -— 



sech 



(46) 

(47) 

(48) 
(49) 



A. Local density of states 

From the quasi-classical Green's function in cq. (|46|) . we can calculate LDOS by eq. (fTTT) as 



for the Majorana bound state \E\ < A B and as 



N(z,E) = N 



\E\ 



\E\ 



\E\ 



■ tan 



VE 2 -A% 2 \y/E*-A% A B ^fW^K\ 



sech ( — - - 



(50) 



(51) 



for the continuum state \E\ > A B . The LDOS at the 
edge z — is shown in Fig. [3{a). This LDOS in the Ma- 
jorana bound state has a linear energy dependence with 
a slope (tt/4)Nq. This linear dependence is also obtained 
by numerical calculation [y, [33|, |34[ • A "Majorana cone" 
formed by the dispersion of Majorana QPs @, Q is the 
cause of this linear dependence. The LDOS in the contin- 
uum state is also deformed to conserve the DOS against 
the appearance of Majorana QPs. 



B. Edge spin current 

The edge mass current is absent in the superfluid 3 Hc 
B-phase. This is clearly seen from eq. because go is 
a real function. We can calculate the spin current using 
cq. (|10l) . The finite spin currents along the edge are 

j x s (z,E) = j s (z,E)(sin9 L x + cos 9 L y), 
f s (z,E) = -j s (z,E)(cose L x-smO L y), (52) 

where x and y are unit vectors. The spin current function 

is{z,E) is 



ith E\E\ l9 / z 

^^--^^seclr 2 ^- 



(53) 



for the bound state \E\ < As and 



1 h E 

3^E) = ---v F N Q — 



^W^A 



E 2 



A B y/E 2 



A 2 



E 2 
3-r^r tan 
A 2 

^B 



A, 



yfW 



A 2 

^B 



sech I — 

SB 



(54) 



for the continuum state \E\ > Ab. This spin current 
function at the edge is shown in Fig. EJb). The en- 
ergy spectrum of the spin current from Majorana QPs 
is quadratic. The edge spin current from the continuum 
state is due to the deformation of the LDOS, which de- 
creases away from Ab. The asymptotic behavior of the 
spin current function in eq. (|54|) at z = is estimated as 



j.(z = 0,E) 



l h 

302 



v F N Q 



Ab 
E 



(55) 



for \E\ 3> Ab, implying that the contribution of QPs 
decreases with the same power law ~ E~ 3 , as in the A- 
phase. This power law behavior of such a decrease means 
that QPs carrying the edge spin current are distributed 
neither in the narrow energy shell around the Fermi level 
nor up to the bottom of the Fermi sea. 

Since QPs fill the energy state up to the Fermi energy 
at a zero temperature, the spin current from the Majo- 
rana bound state is obtained as 



Jf\z) 



f dEj s (z,E) 

J-Ap, 



'-A 

-^ VF N A B sech 2 \J- ) 



and that from the continuum state is 
-A B 



J s cont M 



1 



dEj s (z,E) 
4 \ h 



-v f N Ab sech 



(57) 



24 V 3tt/ 2 Us, 

The total spin current from the Majorana bound state is 



/•oo 



and that from the continuum state is 



(58) 



~3W2^' (59) 



where the quantum of circulation k = hj (2m) emerges 
from the ratio of coefficients between the spin and mass 
currents. Since the edge currents from the Majorana 
bound state and continuum state flow oppositely, sim- 
ilarly to the A-phase, the total spin current by the edge 
state is 



J, = J, 



MJ 



Tcont _ K Ilk 

s ~ 2ir 6 ' 



(60) 




-A B A c 
E 



FIG. 3. (Color online) Energy profiles of LDOS (a) and spin 
current function (b) at z — in the B-phase. At a zero tem- 
perature, QPs fill the colored (shaded) states in (a). The spin 
currents from the bound and continuum states are derived by 
integrating the blue (light gray) and pink (gray) regions in 
(b), respectively. 



The numerical coefficient of nh in the total spin current 
is 2/3 of that in the total mass current because QPs con- 
tributing to the spin current are QPs in two of three spin 
states. Finally, each spin component of the total spin 
current is 

-(sin^a: + cos6lV), 



J ! > 



2tt 6 
k nh 
2^~Q 



(cos9lx — sindiy). 



(61) 



C. Temperature dependence 

The temperature dependence of the total spin current 
J S (T) is calculated using eq. <j8j> as 



UT) 



+2 



Ab^/o 



k nhirkBT 
2?rT A b 

3- 



E 



V u n + A % 



qW 2 7^+Al-A, 



A| 



A 2 



In 



u> n 



(62) 



This temperature dependence is shown in Fig. 0] by the 
open circles compared with the supcrfluid density p° s by 
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FIG. 4. (Color online) Temperature dependence of total spin 
current (open circles) in the B-phase with the superfluid den- 
sity (solid line). 



TABLE I. Comparison between the A- and B-phases for a 
typical d- vector, time reversal symmetry (TRS), a kind of 
edge current, dispersion of Majorana fermions, energy depen- 
dence of DOS in the bound state, energy dependence of edge 
current in the bound state, amplitude of total edge current, 
low-temperature depletion of edge current, and energy states 
where the odd- frequency pair amplitude is related to LDOS. 



A-phase 



B-phase 



d-vector 

TRS 

kind of edge current 

dispersion 

DOS 

edge current 

total edge current 

depletion 

odd-frequency 



broken 

mass current 

Majorana valley 

constant (No/2) 

_E-linear 

\J\ = nh/i 

T 2 -power 

zero energy state 



k x x + k y y + k z z 

unbroken 

spin current 

Majorana cone 

E- linear 

_E-quadratic 

J a \ = (re/27r)(nft/6) 

T 3 -power 

bound states 



the line. These behaviors are clearly different because 
the low-temperature depletion of the total spin current 
is due to the excitations of Majorana QPs. There is no 
other low energy excitation. 

According to cq. (|62l) . we can derive the low- 
temperature behavior of the total spin current as 



UT)—^ 



k nh 
2~^~6~ 



\-C 



irksT 







■nksT 



(63) 



where the coefficient C is fixed within 3/5 < C < 1 by the 
Euler-Maclaurin formula using up to the fourth Bernoulli 
number. This T 3 -power behavior of the depletion comes 
from the excitations of Majorana QPs because j s (z, E) in 
the Majorana bound state has a quadratic energy depen- 
dence [Fig. EXb)] . The observation of the depletion could 
establish the existence of Majorana QPs. Since the deple- 
tion is proportional to T 3 , the low-temperature depletion 
of the spin current in the B-phase is more gradual than 
that of the mass current in the A-phase. 



V. SUMMARY AND DISCUSSION 

We have found analytic solutions of the quasi-classical 
Green's function g in eq. (fTTj) and eqs. (|4"6")) - ([4T)| for the 
superfluid 3 He A- and B-phases with a specular edge, 
respectively These solutions are self-consistent in the 
weak-coupling limit. By using the solutions, the DOS 
including Majorana QPs bound in the edge state has been 
obtained. The mass current in the A-phase and the spin 
current in the B-phase carried by Majorana QPs have 
also been obtained. A comparison between the A- and 
B-phases for the topological features and our main results 
is summarized in Table HI 

In the A-phase confined to a thin slab, the amplitude 
of the total edge mass current is \J\ — nh/A. The angu- 
lar momentum by this edge mass current is L — Nfh/2, 



which corresponds to one of the predictions of the intrin- 
sic angular momentum [29j, |30[. This angular momen- 
tum can be separated into the contributions of Majorana 
QPs, i MJ = Nh, and of the QPs in the continuum state, 
£cont _ _7v^/2. Thus, the total angular momentum 
is given by the cancellation of half of the angular mo- 
mentum from the Majorana bound state by that from 
the continuum state. The cancelled angular momen- 
tum Nh/2 may be related to the DOS of the Majorana 
QPs Nq/2. The reduction in the angular momentum by 
temperature behaves as that in T 2 -power within a low- 
temperature range by the excitations of Majorana QPs. 
Note that, in 3D thick slabs, the angular momentum is 
reduced by the canting of the Z-vector at the side edge [f| . 

In the B-phase, the amplitude of the total edge spin 
current is \J S \ = (K/2n)(nh/6). The numerical coefficient 
of nh for the total spin current is 2/3 of that for the 
total mass current in the A-phase because the QPs in 
two of three spin states contribute to the spin current. 
The cancellation of a portion of the edge current from 
the Majorana bound state by that from the continuum 
state occurs in the B-phase as well as in the A-phase. 
QPs contributing to the edge current are not confined 
in only the vicinity of the Fermi surface, the same as 
those in the A-phase. On the other hand, the reduction 
in the total spin current by temperature behaves as that 
in T 3 -power within a low-temperature range owing to 
the quadratic energy spectrum of the spin current in the 
Majorana bound state. 

Finally, we point out a method of observing Majorana 
QPs. The low-temperature depletion of the edge current 
is due to the excitations of Majorana QPs. Therefore, the 
observation of the depletion of the mass current as that in 
T 2 -power or that of the spin current as that in T 3 -power 
could establish the existence of Majorana QPs. We also 
point out another possibility for experimentally observ- 
ing Majorana QPs by controlling the surface roughness. 
The surface roughness can be controlled by coating the 
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surface with 4 He atoms, which has already been proved 
experimentally |35|]. The edge current from the Majo- 
rana bound state and that from the continuum state can 
be suppressed obeying different rates as a function of 
the surface roughness. Particularly in the B-phase, since 
Majorana QPs are topologically protected, Majorana fea- 
tures are maintained under a certain surface roughness. 
In fact, the linear dispersion of the surface Majorana 
bound state in the supcrfluid 3 Hc B-phase remains intact 
even at the partially rough surface [9|, [K| . The concrete 
calculation for quantitatively estimating this change is a 
future problem. We note that the ed ge m ass currents in 
the 2D A-phase for a diffusive edge |36[ and a retrore- 
flecting edge [18| have been calculated. 



satisfy the gap equation 

A{r,k) = N TTk B T }] (v(k,k')f{r,k',io n ) 



E 



k' 

(Al) 



where lo c is the cutoff energy. The pairing interaction 
V(k,k') = 3gik ■ k' for Cooper pairs with an orbital 
angular momentum 1 = 1, where g\ is the coupling con- 
stant. The anomalous quasi-classical Green's function is 
described in spin space as 



/ = (ia-(Ty) ■ f 



~fx + ify fz 

fz fx + ify 
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Appendix A: Self-Consistency of Solutions 

The self-consistent pair potential A(r, k) and the 
anomalous quasi-classical Green's function f(r,k,io n ) 



1. A-phase 

In the A-phase from the solution of the Riccati ampli- 
tude in eq. (fj~5|) . the spin component of the anomalous 
quasi-classical Green's function f z — 2a/ (1 + ab) is only 
finite in / as 



f z (x,k,u n ) 
x d z (x, k) - 



1 



% + A 2 A sm 2 9 

A 2 , sin 2 9 cos 2 - 



2(oj n + iAa sin 9 sin c 



sech 



u 



(A2) 

The second term of f z contains the odd-frequency pair 
amplitude 0, [13 



ff*(x,k,W n ) 



1 A 2 , sin 2 9 cos 2 . 



A 2 , sin 2 



+ 



1 



uj n + iAa sin # sin u) n — i Aa sin 6 sin 



sech 



u 



(A3) 



The odd-frequency pair amplitude is related to the LDOS, as pointed out by Higashitani et al [38|. In the Majorana 
bound state \E\ < A^ sin#, the real part of the odd-frequency pair amplitude for the energy 



Re [f? F (x,k,E)] =Re [/° F (x, k,cj n )\ iUn ^ E+iri ] 



.r 



- A A Sln flcOS- 0_ ^ E _ ^ gin e sin ^ +6 ( E+ Aa sm Q sin 0)] sech 2 

4 ,/A 2 sin 2 ^-^ 2 ^ A 



(A4) 



and the angle-resolved LDOS 



have the relation 



N(x,k,E) _ n A A sm 2 9cos 2 <f> 
No 2 ^A\sm 2 9-E 2 

x 5(E — Aa sin 9 sin <f>) sech 



N(x,k,E = 0) 

An 



Re 



f OF (x,k,E = Q) 



(A5) 



in our solution with a specular edge. 



(A6) 
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At a zero temperature, the z-componcnt of the d- vector 
is obtained using the gap equation eq. (|A1[) as 



d z (x,k) = giN A A 

■5 



'^ +ln ^) tanh fe 



-ik,. 



6 



In- 



1 



scch 



(A7) 



Since the coupling constant has the relation 



1 5 u c 2u c 

= — h In — — = In — . 

5i A 6 A A Abcs 



(A8) 



where Abcs is the amplitude of the superconducting 
gap for the conventional s-wave superconductivity, the 
^-component of the d-vector is 



d z {x,k) = Aj 



k x tanh ( — — 



+ik v { 1 + 



<7i^o 



scch 



u 



(A9) 



In the weak-coupling limit uj c — > oo or g\ Nq — > 0, 



d z (x,k) = A, 



k r tanh 



t 



■ik, 



(A10) 



in eq. (|13l) is a self-consistent pair potential. 



2. B-phase 

In the B-phase with the solution of the Riccati am- 
plitude in eq. (|43|) . the spin components of the anoma- 
lous quasi-classical Green's function / = M'f'M* are 
described as 



f x =cos((f) + 9 L ) Im 
f y =sin((j) + 6 L ) Im 



2a' 




l+a 
2a' 


'V 


l + a 


b'_ 



fz = Re 



2a' 



1 + a'V 

(All) 



Since 



2a' 



1 + a'b' V^+Al 
x d'(z,8) 



A| cos 2 e 



sech I — - 



2(w n -MA,Bsin6>) ~ \^b 

by an analogy to the A-phase in eq. (|A2|) . we obtain 



(A12) 



f x (z,k,uj n ) 
f y (z,k,u n ) 

fz(z,k,U n ) 



cos(4> + 6l) 

" V"l + A| 

sm((f> + 6 L ) 

1 



V^+A 



Ab sin + i 



A% cos 2 1 



„ A 2 cos 2 6> 
A B sin0 + i— 



1 



1 



w n + iAb sin cj„ — iAb sin ( 
1 1 



Ab cos tanh [ — - ) - 



uj n + i As sin 9 uj n — iAb sin 
A 2 R cos 2 6» / 1 1 



sech 
sech 






(A13) 
(A14) 



scch I 

cj„ + zA_b sin ^ o; n — zA^sin^y \£b 



(A15) 



The second term of f z is the odd-frequency pair amplitude. In the Majorana bound state \E\ < Ab, the real part of 
the odd-frequency pair amplitude for the energy 



Re[f? F (z,k,E)] 



4 y'Ag - E 2 Vsb 



and the angle-resolved LDOS 

N(z, k, E) _ 7T A| cos 2 



_ [ (£ - A S sin0) + S(E + A B sin*)] sech 2 (^- ) 



(A16) 



(A17) 



have the relation 

N(z,k,E) 

Ao 



Re 



r"{z,k,E) 



in our solution with a specular edge. The odd-frequency 

pair amplitude is related to the LDOS not only in the 

(A18) zero-energy state as in the A-phase but also in the Ma- 
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jorana bound state. 

At a zero temperature, each component of the d- vector 
is obtained using the gap equation eq. (|A1[) as 



d x {z,k) =giN A B (k x cos6 L - k y smO L ) 

2w c 1 2 ( z 

In — - + - sech — 
A B 6 \£b 

d y (z,k) =giN A B (k x sm6L + ky cos 6£) 

2uj r 1 



In- 



- sech 2 I — - 
6 Ub 



d z (z,k) =g 1 N Q A B k z \n-— tanh — - 



Since the coupling constant has the relation 



1 2uj c 2lu c 

— ^- = In -7 — = m . 

giN A B Abcs 



(A19) 

(A20) 
(A21) 



(A22) 



each component of the d- vector is 

d x (z,k) =A B (k x cos9 L -k y sm0 L ) 





X 


", , 31^0 ,2 ( * Y 


dy(z,k) 


=A B (k x sm0L + k y cos9l) 




X 


) +! ¥ -*(£)] 



d z (z,k) =A B k z t&nh — 
,£b 



(A23) 

(A24) 
(A25) 



Bumps at the edge in the k x - and /^-components are 
smaller than that in the fey-component for the A-phasc 
according to the coefficients of the squared hyperbolic 
secant function. In the weak-coupling limit uj c — > oo or 
3i^o -> 0, 



A(z,fe) = A i 



iin0e- l ^ +e ^ cos0tanh(z/6?) 



cos0tanh(z/£s) sin#e 



i{4>+e L ) 



(A26) 



in eq. (|37| is the self-consistent pair potential. 
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